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Abstract 

We derive the exact formula for thermal-equilibrium spacing distribution of one-dimensional particle gas 
with repulsive potential depending on the distance r between the neighboring particles. We are focused on 
the power-law potentials r~ a for a > 0, being motivated by the actual traffic research. 



PACS numbers: 05.70-a, 05.20-y 

Investigation of one-dimensional particle ensembles 
seems actually to be very useful for understanding of 
the complex system called vehicular traffic. Beside 
the popular cellular automata, in the recent time new 
trend appears in the traffic modelling. Application of 
the equilibrium statistical-physics to non-equilibrium 
traffic systems has been successfully demonstrated for 
example in the Ref.Q]. In this continuing work we 
analyze a certain family of the particle gases exposed 
to the heat bath with the temperature T > 0, being 
primarily aimed at the equilibrium distribution for 
particle distance (spacing distribution) and particle 
velocity (velocity distribution). Subsequently, such a 
detailed description can be used for comparing with 
the microscopic structure of the realistic traffic flows 
(in Ref.[U). 

Consider N identical particles on the circle of the 
circumference L = N. Let X4 (i = 1,2,..., TV) rep- 
resents the circular position of the i— th particle. Put 
xn+i = x\ + 2ir, for convenience. Now we introduce 
the short-ranged potential energy 



JV 
i=l 



v(n), 



tween the neighboring particles only. Thus, the hamil- 
tonian of the described ensemble reads as 



N 



N 



n=lY,(v t -v) 2 + CY / V(r 1 ), 

i=i i=i 

with the i— th particle velocity Vi and the positive con- 
stant C. Note that v represents the mean velocity in 
the ensemble. Then, the appropriate partition func- 
tion 1 



dr\dv\ . . . dr^dvN 



(1) 

leads us to the simple assertion that velocity v of par- 
ticles is Gaussian distributed, i.e. 



P(v) 



1 



e 2kT 



is the corresponding probability density. 

Of larger interest, however, is the spacing distri- 
bution Pp(r). In order to calculate the exact form of 
Pp (r) , one can restrict the partition function to the 
reduced form 



where V(ri) corresponds to the repulsive two-body po- Zn(L) = 
tential depending on the distance r,; = |x{+i— Xi\j- be- 
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1 k represents the Boltzmann factor 



dr 



N 



1 



where (3 — is so-called inverse temperature of the 
heat bath. Denoting f(r) — e~P v v> the previous ex- 
pression changes to 

Z N (L) = [ 51 L Tr,] TT/( r <) dr i ■ ■ ■ dr N- 
^ N \ ~{ J t i 

Applying the Laplace transformation (see the Ref.[2j 
for details) one can obtain 

/>OG 

g N (p)= / Z N (L) e- pL dL = 
Jo 



f(r)e-P r dr\ = [g(p)} 



N 



iN 



Then the partition function (in the large N limit) can 
be computed with the help of Laplace inversion 



Z N {L) 



1 

2~7ri 



B+ioo 



g N (p) e Lp dp. 



Let us to proceed to the special variants of the gas 
studied. Firstly, we draw our attention to the Coulomb 
gas with the logarithmic potential 

V{r) := -ln(r). 

Such a gas is frequently used in the many branches of 
physics (including the traffic research in Ref. 0]) and 
the corresponding spacing distribution reads as (see 
Ref. [2]) 



r(/? + i) 



B — ioo 



where r(£) is gamma function. Of larger physical inter- 
est, as demonstrated in Ref. [I] and PJ, seems actually 
to be the potential 

V a (r) :=r~ a , 

for a > 0. The aim of the following computational 
procedure is to normalize the distribution 



Its value is well estimated by the approximation in the 
saddle point B which is determined using the equation 

1 dg L 
g(B)d P [ ' N' 



Thus, 



Z N (L) « [g(B)]" e 



N LB 



(2) 



Hence the probability density for spacing r\ between 
the particles j}2 and jjl can be then reduced to the form 

P(ri) = — — fin). 

Z N {L) 

Supposing N 1 and using equation <|2J we obtain 

P(n) = 4r/(n) e- Br \ 

which leads to the distribution function for spacing r 
between arbitrary couple of neighboring particles 



P p (r) = A e-^W e- Br . 



(3) 



Note that constant A assures the normalization 
J °° Pp(r) dr = 1. Furthermore, returning to the orig- 
inal choice L = N, the mean spacing is 



(/■) - z j rPp(r) dr = 1. 
'o 



(4) 



Two above conditions can be understood as equation 
system for unknown normalization constants A, B. 



Pp(r) = Ae 



-Br 



(5) 



Consider now the favorable choice a = 1, for which 
the normalization integrals are exactly expressed as 



J e-^e- Br dr = 2^ -| £i (2y/0) 



re 



-Br 



dr = 2^K 2 [2^pB 



(6) 



where JC\ is the Mac-Donald's function (modified 
Bessel's function of the second kind) of order A, having 
for A = 1 and A = 2 an approximate expression 



KM = ^ V (V 1/2 + ^ V 3/2 + O (y- 5 ^ . 

Applying the equations (HI to the normalization inte- 
grals one can determine the exact values of the con- 
stants A and B. Both of them can be, moreover, very 
well estimated by the large (3 approximations 



B 



and 



900 (3 



(3 - 16/5 + ^9 + 864/3 + 256/?< 



WKi (^2/3(2/3 + 3) 



2 



Finally, we investigate the distribution © for gen- 
eral a > 0. Although in this case the normalization 
integrals are not trivially solvable, the scaling (@J leads 
us to the simple formula 

CI 

B « a/3 + 1 + - (0->oo). (7) 

The rough estimation r~ a w 1 — a + ar^ 1 provides the 
asymptotical formula for normalization constant A : 
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Figure 1. Relative deviation in the approximate value of the 

normalization constant B as function of parameter j3. 
We display the deviation 1101 between the numerical value B ex 
and the value B es t, obtained from the large j3 approximation 

J7J. The plus signs, squares, triangles, circles and stars 
correspond to the parameters a = 1,2,3,4,5, respectively. The 
tails of the curves are magnified on the insert. 

For practical applications it seems to be useful 
to detect the critical inverse temperature f3 cr u under 
which the relative deviation between the exact (ex) 
and estimated (est) values of the constant A (or B) 

_ \log(A ex) -log(A est )\ 



S B ■= (10) 

are larger then the fixed permissible deviation 6. For 
these purposes we plot the functional dependence 5b = 
Sb(P) and $]o & (A) = <$iog(A)(/3) in the Figure 1 and 
Figure 2, respectively. We note that the exact values 
A ex , B ex were determined with the help of the numer- 
ical computations. 
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Figure 2. Relative deviation in the approximate value of the 

normalization constant A as function of parameter f). 
Plotted is the deviation J^J between the numerically computed 
value A ex and the estimated value J8j. The symbols used here 
are consistent with the symbols in the Fig. 1. 

To conclude, we have found the exact form of 
the thermal-equilibrium spacing distribution for one- 
dimensional gas which neighboring particles are re- 
pulsed by the two-body potential V = r~ a , where r 
is their mutual distance. The values of two normaliza- 
tion constants were successfully estimated by the sim- 
ple approximation. The determined distribution has 
been and will be applied in the traffic research where 
compared to the distance clearance distribution of the 
freeway traffic samples. 
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